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/  FOREWORD 

£ 

When  the  GEOSAT  satellite  was  launched  on  March  12, 1984,  it  carried  onboard  a  radio  altimeter  that 
continuously  generated  measurements.  The  Naval  Surface  Weapons  Center  (NSWC)  was  given 
responsibility  for  the  reduction  of  this  data  to  best-estimated  along-track  geoid  heights  and  vertical 
deflections.  7/  ■> ,  r  ■,  s  ,  ■  c-  -Vv  U  . -  J-> 

These  geodetic  results  are  currently  distributed  to  the  Defense  Mapping  Agency  Aerospace  Center 
(DMAAC),  St.  Louis,  Missouri;  the  Naval  Oceanographic  Office  (NAVOCEANO),  Bay  St.  Louis,  Mississippi; 
and  other  specified  users  for  geodetic  and  oceanographic  analyses. 

This  technical  report  was  reviewed  and  approved  by  Ralph  L.  Kulp,  Jr.,  Head,  Space  and  Ocean 
Geodesy  Branch,  and  by  Carlton  W.  Duke  Jr.,  Head,  Space  and  Surface  Systems  Division. 
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1.0  INTRODUCTION 


The  GEOSAT  mission  is  comprised  of  a  primary  and  secondary  mission.  The  primary  mission 
meets  the  geodesy  requirement  to  provide  the  Department  of  the  Navy  (DON)  with  a  global  data  base 
of  10  cm  precision  radar  altimeter  measurements.  This  data  will  allow  for  improvements  in  the 
gravitational  models  required  by  advanced  Submarine  Launched  Ballistic  Missile  (SLBM)  systems. 
The  secondary  mission  meets  the  oceanography  requirement  to  provide  the  DON  with  altimeter  data 
for  the  GEOSAT  Ocean  Applications  Program  (GOAP).  Several  agencies  are  responsible  for  the 
GEOSAT  mission. 

The  Naval  Surface  Weapons  Center  (NSWC)  has  as  a  responsibility  the  development  of  an 
analysis  system  for  the  reduction  of  radar  altimeter  measurements  to  the  best-estimated  along-track 
geoid  heights  and  vertical  deflections.  These  geodetic  products  are  to  be  distributed  to  the  Defense 
Mapping  Agency  Aerospace  Center  (DMAAC),  St.  Louis,  Missouri,  the  Naval  Oceanographic  Office 
(NAVOCEANO),  Bay  St.  Louis,  Mississippi,  and  other  specified  users. 


2.0  GEOSAT  PROJECT  DATA  FLOW 


Figure  1  illustrates  the  flow  of  data  through  the  GEOSAT  project ,  showing  the  source,  data  file 
type,  and  user(s)  of  the  major  data  files.  The  Johns  Hopkins  University  Applied  Physics  Laboratory 
(JHU/APL)  is  responsible  for  implementation  of  the  altimeter,  spacecraft,  control  center  and  ground 
station.  Functions  that  are  performed  at  the  Control  Center  include  health  monitoring,  command 
generation,  and  data  management.  Data  collected  by  the  altimeter  and  other  spacecraft  functions  are 
stored  on  tape  recorders  and  dumped  to  the  APL  ground  station  every  12  hr  on  the  S-band  downlink  to 
the  station’s  60-ft  dish.  In  addition  to  controlling  the  satellite,  APL  archives  the  raw  data,  uses 
station  alerts  from  Defense  Mapping  Agency  Hydrographic/Topo graphic  Center  (DMAHTC)  and 
transmits  preprocessed  data  to  NSWC,  the  Naval  Research  Laboratory  (NRL),  and  the  Naval 
Research  and  Development  Activity  (NORDA).  DMAHTC  employs  the  standard  DMA  Doppler- 
combined  Tranet/Magnavox  1502  D.S.  tracking  station  network  in  collecting  Doppler  Data. 


2.1  SENSOR  DATA  RECORD 


The  Sensor  Data  Record1  (SDR)  is  the  data  interface  between  APL  and  NSWC.  It  consists  of 
radar  altimeter  data  with  instrument  and  spacecraft  system  errors  identified  and  quantified.  NSWC 
began  receiving  SDR  tapes  via  messenger  on  a  regular  schedule  after  the  satellite  was  stabilized  and 
the  APL  data  processing  procedures  became  operational.  The  SDR  is  the  primary  input  to  the  radar 
altimeter  data  processing  system  at  NSWC. 
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□  -  SOURCES  O-  PROCESS  =■  DATA  FILES  — ►  -  DATA  FLOW 


FIGURE  1.  GEOSAT  PROJECT  DATA  FLOW  DIAGRAM 
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i  2.2  FILTERED  GEOPHYSICAL  DATA 


Filtered  Geophysical  Data  (FGD)  files  produced  at  NSWC  are  mailed  to  DMAAC,  NAVOCEANO, 
and  other  specified  users.  The  FGD  consists  of  best-estimated  along-track  geoid  heights,  vertical 
deflections,  sea  surface  heights,  environmental  corrections,  statistics,  and  data  quality  indicators. 
The  format  is  given  in  Section  10.3. 


2.3  WAVEFORM  DATA  RECORD 


The  Waveform  Data  Record1  (WDR)  tape  is  the  data  interface  between  APL  and  NRL.  It  contains 
waveform  samples  taken  at  10-points/sec  (pps)  intervals.  NRL  applies  this  information  in  conducting 
radar  altimeter  validation  studies. 


2.4  INTERMEDIATE  GEOPHYSICAL  DATA  RECORD 


In  addition  to  the  WDR  from  APL,  NSWC  sends  NRL  Intermediate  Geophysical  Data  Record 
(IGDR)  tapes  for  ocean  and  land/ice  data  produced  at  the  10-pps  rate.  These  tapes,  detailed  in  Section 
10.2,  are  output  from  the  preprocessor  program.  In  addition  to  the  altimeter  measurements,  the 
IGDR  includes  the  precise  orbit  and  environmental  corrections. 


2.5  NAVAL  OCEAN  RESEARCH  AND  DEVELOPMENT 
ACTIVITY  DATA  RECORD 


To  satisfy  the  GOAP  requirement,  Naval  Ocean  Research  and  Development  Activity  (NORDA)’s 
Data  Record  (NDR)  is  transmitted  directly  from  APL  via  data  line  to  NORDA,  where  it  is  employed 
for  oceanographic  studies.  At  a  later  point  in  the  mission,  the  NDR  may  be  transmitted  from  APL 
directly  to  Fleet  Numerical  Oceanography  Center  (FNOC),  Monterey,  California. 


2.6  UNCLASSIFIED  DATA  RECORD 


Possible  distribution  of  an  Unclassified  Data  Record  (UDR)  for  civilian  use  is  currently  being 
discussed.  The  specific  content  and  methods  of  creating  the  UDR  have  not  yet  been  defined 


2.7  DOPPLER  DATA  RECORD 


In  addition  to  its  altimeter,  the  satellite  carries  a  Doppler  beacon  that  transmits  data  to  a 
worldwide  network  of  Doppler  stations  DM  AHTC  collects  and  preprocesses  (cleans)  the  data,  and 
then  sends  it  to  NSWC  for  use  in  generating  precise  orbits,  one  of  the  major  inputs  into  the  NSWC 
altimeter  data  processing  system 
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3.0  RADAR  ALTIMETER  DATA  PROCESSING  SYSTEM 


3.1  GENERAL  DESCRIPTION 

The  radar  altimeter  data  processing  system  reduces  the  radar  altimeter  measurements  to  best- 
estimated  along-track  geoid  heights  and  vertical  deflections.  This  is  accomplished  by 

computing  environmental  corrections, 

preprocessing  altimeter  measurements, 

generating  sea  surface  heights, 

computing  geodetic  position  and  velocity  along  track, 
smoothing  the  raw  geoid  heights, 
archiving,  and 

distributing  the  final  geodetic  products  to  specified  DMA  users. 

In  generating  the  final  geodetic  products,  the  processing  system  uses  several  input  data  bases: 
world  map  coastline  coordinates, 
environmental  data, 
precise  orbit, 

tide  amplitudes  and  phases, 

the  satellite  revolution  (rev)  epoch  table,  and 

the  SDR 

Figure  1  shows  the  data  flow  for  the  system  and  user  data  interfaces.  The  environmental 
corrections  generator,  preprocessor,  segmenter,  and  smoother  comprise  the  major  software 
components. 


3.2  DATA  FLOW  OF  PROCESSING  SYSTEM 


The  data  flow  diagram  for  the  NSWC  altimeter  data  processing  system,  Figure  2,  shows  major 
components  of  the  NSWC  radar  altimeter  data  processing  system  and  data  interfaces  The  SDR, 
world  map  coastline  coordinates,  precise  orbit,  tidal  constants,  satellite  rev  epoch  table,  sunspot 
number,  and  climatology  data  comprise  the  major  inputs  to  the  NSWC  radar  altimeter  data 
processing  system  Output  files  produced  by  the  system  are  the  ocean  IGDR,  land/ice  IGDR,  FGD,  and 
plots 
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The  classified  SDRs  generated  at  the  JHL'/APL  are  written  to  tape  and  transported  by  messenger 
to  NSWC,  White  Oak,  Maryland,  and  then  to  NSWC,  Dahlgren,  Virginia.  The  SDR  tapes  are  9-track 
(1600  bpi)  and  data  are  written  in  ASCII  format.  The  newly  installed  CDC/875  computer  is  used  in 
the  reduction  of  the  altimeter  measurements  to  along-track  geoid  heights  and  vertical  deflections. 
Remote  terminals  are  available  for  accessing  the  CDC/875,  a  secure  computer. 


FIGURE  2.  NSWC  GEOSAT  ALTIMETER  DATA  PROCESSING  SYSTEM 

DATA  FLOW 


Three  world  map  data  bases  are  available  for  separating  the  land  and  ocean  data:  Hershey 
WorldO  (8000  data  points),  CIA  Worldl  (60,000  data  points);  and  NAVOCEANO  World2  (1,700,000 
data  points).  World2  contains  the  detail  required  by  the  GEOSAT  project  in  the  separation  of  land 
and  ocean  altimeter  measurements.  A  masking  grid  was  derived  from  the  World2  data  base  It 
consists  of  land/ocean  latitude  and  longitude  boundaries  for  determining  whether  a  specific 
measurement  is  over  land  or  ocean.  The  land/ice  data  tape  also  contains  data  that  has  failed  ocean 
data  tests  criteria.  These  data  are  referred  to  as  "other”. 


The  precise  orbits  (long  arcs)  are  computed  in  2-day  time  spans  using  observations  from  the 
combined  Tranet/Magnavox  1502  D  S.  Doppler  tracking  network.  The  precise  orbits  are  smoothed  at 
two-day  crossovers  and  stored  on  the  CDC/875  computer  as  permanent  files.  The  CELEST  computer 
program  for  orbit  computations  and  file  formats  are  described  by  O’Toole2.  The  orbit  file  names  have 
a  fixed  format:  N17YDDD,  where  N  is  either  "L”  or  "S",  designated  long  arcs  or  long  arcs  smoothed  at 
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two-day  crossovers  respectively.  The  GEOSAT  satellite  identification  number  is  17;  Y  is  last  digit  of 
year,  and  DDD  denotes  a  3-digit  day  number  of  first  day  of  two-day  fit. 

The  Schwiderski  ocean  tide  model3  requires  harmonic  tidal  constants  (amplitudes  and  phases) 
along  a  standard  GEOSAT  satellite  track.  The  amplitudes  and  phases  of  the  leading  harmonic  ocean 
tidal  modes  (M2,  S2,  N2,  K2,  Kl,  01,  PI,  and  Ql  with  solid  earth  tide  and  loading  effects  included) 
are  stored  as  permanent  files  and  used  to  determine  the  tidal  effects  on  the  altimeter  measurements. 
Other  input  to  the  model  includes  year,  day,  sec,  time  interval,  total  number  of  points,  latitude, 
longitude,  longitude  of  ascending  node,  and  the  equator-crossing  time. 

The  satellite  rev  epoch  table  consists  of  year,  day,  rev  number,  longitude  of  the  ascending  node, 
time  of  the  ascending  node,  an  average  (over  2-day)  period,  and  an  average  (over  2-day)  shift  of  the 
longitude  of  the  node.  These  values  are  computed  using  precise  orbits,  and  are  required  because  of 
the  noncircularity  of  the  GEOSAT  orbit.  The  table  entries  are  given  for  every  100th  rev  number. 

Temperature,  atmospheric  pressure,  and  partial  pressure  of  water  vapor  are  the  climatology  data 
fields  transmitted  daily  from  FNOC  via  the  Naval  Oceanographic  Data  Distribution  System  to  the 
NSWC  Mohawk  data  communications  computer.  These  parameters  are  acquired  in  two  12-hr 
periods/day  and  archived  for  use  in  the  computation  of  barotropic,  wet  tropospheric,  and  dry 
tropospheric  conditions.4- 5 

Output  files  from  the  data  analysis  and  reduction  system  include  land/ice  and  ocean  IGDRs,  FGD, 
and  plots.  The  land/ice  data  IGDRs  are  archived,  sent  to  NRL  in  Washington,  D  C.,  and  applied  in- 
house  for  topographic  analyses.  The  ocean  IGDRs  at  2-pps  are  archived  and  used  as  the  main  input  to 
the  segmenter. 

System  2000  was  initially  selected  as  the  data  base  management  system  for  maintaining  an 
accurate  record  of  the  GEOSAT  data  files  and  related  information.  However,  after  this  selection  was 
made,  a  Relational  Information  Management  (RIM)  data  base  management  system  became  available, 
a  system  that  better  satisfied  the  project’s  needs.  Input  to  RIM  consists  of  a  Revs  Table  and  a  Days 
Table.  The  Revs  Table  is  automatically  updated  with  current  information  from  the  FGD  file,  while 
the  Days  Table  is  updated  by  manually  created  files.  Section  6.0  provides  a  description  of  the  data 
base  management  system  elements. 

The  IGDR,  the  output  from  the  preprocessor  program,  becomes  input  for  the  segmenter  program. 
The  file  name  formula  is  XYDDDBP,  where  X  is  an  O  or  L  that  represents  ocean  data  or  land/ice  data 
respectively;  Y  is  the  year,  DDD  is  the  day  number;  B  indicates  binary  information;  and  P  is  a  1  or  2 
representing  10-  or  2-data  pps  respectively. 

Continuous  graphic  output  makes  possible  the  validation  and  verification  of  geodetic  products. 
Optional  on-line  plotting  is  available  for  displaying  parameter- versus-time  relationships.  A 
parameter  may  be  selected  from  several  contained  in  the  geophysical  data,  such  as  significant  wave 
height  (SWH),  automatic  gain  control  (AGC),  altitude,  or  sea  surface  height.  During  production 
processing,  geoid  heights  and  vertical  deflections- versus-time  are  plotted  for  each  rev.  Plots  may  also 
be  viewed  on  the  CRT  via  Tektronix  equipment.  Other  graphic  output  can  be  accessed  from  a  world 
map  computer  program  called  GEOMAP  that  displays  geophysical  data  on  a  world  map  background.4 
A  preprocessor  program  generates  the  GEOSAT  data  file  formats  that  are  compatible  with  those 
required  by  GEOMAP.  The  GEOMAP  program  was  developed  for  the  SEASAT  project  by  Sperry 
Univac  at  Dahlgren,  Virginia. 

The  FGD  data  base  is  organized  sequentially  by  rev.  In  addition  to  the  best-estimated  along- 
track  geoid  heights  and  vertical  deflections,  it  includes  sea  surface  height,  environmental  and  system 
corrections,  quality  flags,  and  statistics.  The  FGD  file,  in  binary  format,  is  converted  to  ASCII  format 
before  its  distribution  to  DMAAC,  N  AVOCEANO,  and  other  specified  users  Since  requests  for 
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smaller  quantities  of  data  are  identified  by  rev,  area,  or  time,  an  auxiliary  computer  program  was 
written  to  extract  data  from  the  data  base  in  fulfillment  of  these  needs. 


3.3  MAJOR  SYSTEM  COMPONENTS 

3.31  Environmental  Corrector 


The  altimeter  measurement  is  determined  by  calculating  the  time  it  takes  to  transmit  the  radar 
signal  from  the  spacecraft  to  the  ocean  surface  and  to  reflect  back  to  the  spacecraft.  The  radar  signal 
will  encounter  environmental  effects  during  the  transit.  The  measurement  also  contains  errors 
caused  by  instrument  and  spacecraft  biases.  Corrections  for  these  errors  are  described  in  Section 
5.32. 

Corrections  for  environmental  effects  are  computed  by  the  Environmental  Corrector  (EN  VCOR), 
taking  as  input  the  precise  orbit,  sunspot  number  table,  rev  epoch  table,  climatology  fields,  World2 
data  base,  tide  amplitudes  and  phases,  and  the  SDR  header.  Output,  written  to  a  temporary  file, 
consists  of  the  ionospheric,  wet  tropospheric,  dry  tropospheric,  barotropic,  and  tide  corrections  for 
data  points  over  the  oceans.  The  corrections  are  computed  for  ocean  data  that  exist  in  the  time  span 
specified  in  the  SDR  header  record.  Ocean  data  points  are  determined  by  a  land  check  routine  that 
uses  a  masking  grid  generated  from  the  World2  data  base  coastline  coordinates. 

Subsatellite  position  at  1-sec  intervals  are  obtained  by  an  8-point  Lagrangian  interpolation  in  the 
precise  orbit  (given  at  60-sec  intervals).  The  rev  number,  time  of  the  ascending  node,  and  longitude  of 
the  ascending  node  are  computed  for  use  in  the  tide  correction  procedure 

Wet  tropospheric,  dry  tropospheric,  and  barotropic  corrections  require  surface  temperature, 
atmospheric  pressure,  and  partial  pressure  of  water  vapor  obtained  from  the  FNOC  climatology  data 
base.  The  environmental  corrections  file,  output  from  EN  VCOR,  is  input  to  the  preprocessor 
program. 

3.32  Preprocessor 

The  altimeter  measurements  as  acquired  from  APL  consist  of  data  taken  over  the  earth’s  surface, 
including  land,  water,  and  ice.  Measurements  may  also  include  spurious  or  wild  data  points.  Before 
the  measurements  can  be  put  into  the  form  required  by  the  Kalman  smoother  for  reduction  to  best- 
estimated  along-track  geoid  heights  and  vertical  deflections  (the  final  NSWC  geodetic  product),  the 
preprocessor  component  must  perform  a  number  of  operations.  These  major  operations  performed  are 
computation  of  the  subsatellite  position  and  velocity  at  measurement  time,  tagging  (and  separation) 
of  land/ice  and  ocean  data,  correlation  of  environmental  corrections,  elimination  of  spurious  data 
points,  aggregation  of  data  points,  and  separation  of  data  into  revs.  Input  to  the  preprocessor  consists 
of  the  SDR,  Doppler  precise  orbit,  environmental  corrections,  rev  epoch  table,  World2  data  base 
coastline  coordinates,  and  tide  amplitudes  and  phases. 

Subsatellite  position  (geodetic  latitude  and  longitude)  and  velocity  at  the  altimeter  measurement 
time  are  computed  by  employing  an  8-point  Lagrangian  interpolation  routine.  The  masking  grid 
referred  to  in  ENVCOR  is  also  employed  by  the  preprocessor  in  the  separation  of  ocean  and  land/ice 
data 

A  10-point  sliding  linear  least-squares  fitting  procedure  eliminates  spurious  data  points  (see 
Section  5.14)  The  SWH  and  AGC  are  similarly  processed  The  altimeter  measurements  given  at  10- 
pps  are  averaged  to  produce  data  at  2  pps  The  number  of  points  averaged  is  provided  as  an  input 
parameter 
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Each  SDR  tape  contains  24-hr  altimeter  data  that  the  preprocessor  separates  into  revs  and  time- 
continuous  segments.  After  this  point  in  the  system,  data  are  processed  according  to  rev,  with  the  rev 
number  as  a  primary  identifier. 

3.33  Sea  Surface  Height  Generator  (Segmenter) 

The  segmenter  uses  the  Doppler  precise  orbit,  altimeter  range  (height)  and  subsatellite  position 
to  compute  the  separation  of  the  reference  ellipsoid  and  the  geoid.  This  separation  is  called  the  sea 
surface  height.  Raymond  Manrique  of  NSWC  Dahlgren's  Space  and  Ocean  Geodesy  Branch  prepared 
the  algorithm  for  computing  subsatellite  position  and  sea  surface  height.6  The  sea  surface  heights 
are  corrected  for  environmental  and  systematic  effects  as  shown  in  Section  5.3.  Output  from  the 
segmenter,  the  Unfiltered  Geophysical  Data  (UGD)  file,  is  input  to  the  smoother  section. 

3.34  Kalman  Smoother  for  Sea  Surface  Heights  (Smoother) 

The  sea  surface  (raw  geoid)  heights  are  reduced  to  best-estimated  along-track  geoid  heights  and 
vertical  deflections  by  an  adaptive  nonsteady-state  Kalman  smoother  This  smoother  is  based  upon 
the  third-order  Markov  process7,  adapted  for  the  altimeter  by  Dr.  C.J.  Cohen  and  described  in  West.6 
An  adaptive  algorithm  by  Ugincius8  provides  improvements  to  the  set  of  a  priori  smoother 
parameters;  autocorrection  distance,  geoid  height  variance,  and  the  noise  standard  deviation. 

The  UGD  file,  input  to  the  smoother,  consists  of  ocean  data  for  one  satellite  rev,  divided  into  time- 
continuous  segments.  Each  time-continuous  segment  is  comprised  of  records,  one  for  each  sea  surface 
height  value.  There  is  a  constant  time  interval  between  records.  An  upper  limit  of  101  min  (period  of 
satellite)  or  approximately  12,000  records  is  used  to  determine  computer  storage  for  one  time- 
continuous  rev  segment.  The  UGD  file  contains  data  at  the  2-pps  rate  during  the  production  mode  of 
operation. 

The  time-continuous  data  segments  are  divided  into  sections  of  approximately  150  sec  (or  1 ,000 
km);  for  each  section  a  cubic  fit,  fit  derivative,  filter  parameters  and  differences  between  the  raw 
geoid  heights  and  the  corresponding  fit  values  are  also  computed.  These  differences,  fit  values,  and 
derivatives  are  saved.  Data  sections  that  are  approximately  1,000  km  in  length  and  having 
approximately  the  same  autocorrelation  distance,  are  grouped  together.  These  are  called  clusters. 

These  data  clusters  are  input  to  the  Kalman  smoothing  procedure,  that  filters  both  forward  and 
backward.  For  each  raw  geoid  height  value  in  a  cluster,  a  smoothed  (or  best-estimated)  geoid  height 
value  and  a  vertical  deflection  value  are  computed.  Any  dubbed-in  data  is  de-weighted  by  setting  the 
reciprocal  of  the  observation  error  equal  to  zero. 

The  FGD  file  is  the  primary  output  from  the  smoother  program  and  includes  the  best-estimated 
along-track  geoid  heights  and  vertical  deflections.  A  detailed  description  of  the  FGD  file  is  given  in 
Section  10.3. 

3.35  List  and  Plot  Programs 

The  list  and  plot  programs  are  controlled  by  LPGEO,  a  procedure  that  produces  a  menu.  From 
this  menu,  list  and  plot  options  can  be  selected  The  list  programs  have  several  options  available  for 
obtaining  printed  output  from  SDR,  IGDR,  and  FGD  files 

SDR  List  Programs  Options 

1.  List  the  entire  SDR  file  or  a  subset  of  the  file 

2.  Send  output  to  the  printer,  laser,  or  microfiche 
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3.  Print  subset  of  file  using  a  start  and  stop  frame  count. 

4.  List  ail  parameters  in  a  data  record  or  a  predefined 
subset  of  the  parameters  in  a  data  record. 

Potions  for  the  IGDR  List  Programs  Printout 

1 .  Print  output  data  on  the  printer  or  laser. 

2.  Print  all  parameters  in  the  data  record  or  a  predefined 
subset  of  the  parameters. 

3.  Print  the  entire  data  file  or  data  identified  by 
rev  number,  time,  or  area 

The  FGD  list  programs  have  as  input  "packed”  and  "unpacked”  data  files  .  Unpacked  files  have 
one  parameter  per  computer  word.  In  order  to  conserve  computer  storage,  these  files  are  used  to 
generated  "packed”  files  that  contain,  where  possible,  several  parameters  per  computer  word. 

The  programs  give  printed  output  from  the  FGD  files. 

Printed  Options  from  the  FGD  Files 

1 .  Print  all  parameters  in  a  record  or  a  predefined 
subset  of  parameters. 

2.  Print  entire  rev  of  data,  or  a  specified  time  span,  or  data  through  a 
specified  area. 

3.  Print  every  nth  data  record. 

4.  Print  output  on  printer,  laser,  or  CRT  (remote  terminal  screen). 

The  SDR  plotting  programs  use  the  major  and  minor  frame  counts  (computes  corresponding 
times)  as  points  of  reference  for  plotted  output.  Plots  may  be  generated  for  all  or  any  of  the  following: 
SWH  ,  altimeter  height,  and  AGC  . 


SDR  Plotting  Program  Options 

1.  Plot  entire  SDR  file  or  a  subset  of  the  file  If  subset,  specify  "start 
major  and  minor  frame  counts”  and  "stop  major  and  minor  frame  counts” 

2.  Plot  all  data  on  one  plot  or  a  specified  data  span  (given  in  sec)  per  plot 

3.  Plot  1, 2,  5,  or  10-pps  data  files. 

4.  Plot  altimeter  height,  SWH,  and/or  AGC. 

The  IGDR  plotting  programs  provide  optional  plots  for  1 1  parameters  They  are  altimeter  height, 
reference  ellipsoid  height,  SWH,  AGC,  sea  surface  height,  ionosphoric  correction,  dry  tropospheric 
correction,  wet  tropospheric  correction,  barotropic  correction,  tide  correction,  and  altimeter  height 
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and  reference  ellipsoid  height  on  same  plot.  The  entire  binary  IGDR  file,  or  one  rev  from  the  binary 
IGDR  file  may  be  plotted.  The  details  of  these  two  options  are  as  follows. 

IGDR  Plotting  Program  Options 

1.  Plot  1  rev  of  binary  IGDR  File 

0  Plot  whole  rev  or  1  time  interval  of  rev  (if  interval,  give  start  and 
end  times). 

0  Select  parameters  to  be  plotted. 

0  Plot  whole  rev  on  one  plot  or  a  specified  number  of  sec  per  plot. 

0  Plot  1  or  2  pps  from  a  2  pps  IGDR  file 
0  Plot  1, 2,  5,  or  10  pps  from  a  10  pps  IGDR  file. 

2.  Plot  entire  binary  IGDR  file  (batch)  laser  plotter. 

0  Select  parameter(s)  to  be  plotted.  The  list  of  1 1  parameters  are 
given  above. 

0  Plot  1  rev/frame  or  a  specified  number  of  sec  per  frame 
0  Plot  from  a  2  or  10  pps  IGDR  file 
o  Plot  1  pps,  or  2  pps  from  a  2  pps  IGDR  file 
0  Plot  1,  2,  5,  or  10  pps  from  a  10  pps  IGDR  file 

3.  Plot  1  FGD  file  (batch/interactive) 

0  Output  plots  to  laser  or  CRT 

0  Plot  entire  rev  on  1  frame  or  a  specified  number  of  sec/frame 

0  Give  start  and  end  sec  if  a  time  interval  has  been  selected 

0  Select  plot  parameter(s),  geoid  height;  vertical  deflection, geoid 
height  and  corrected  sea  surface  height;  and  geoid  height  and  sea 
surface 

4.  Plot  many  FGD  file  (batch/interactive) 

0  Select  output  devices  (fiche,  laser,  or  CRT) 

0  Plot  entire  rev  on  1  frame  or  a  specified  time  span  on  each  frame 

°  Select  parameter(s).  geoid  height,  vertical  deflection,  geoid  height 
and  corrected  sea  surface  height,  or  geoid  height  and  uncorrected  sea 
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4.0  PRECISE  DOPPLER  ORBIT 


NSWC’s  GEOSAT  effort  is  divided  into  three  major  tasks  that  are  performed  in  parallel:  the 
precise  Doppler  orbit  computation;  altimeter  data  analysis  and  reduction;  and  calibration  and 
validation  of  the  geodetic  products.  This  section  describes  briefly  the  precise  Doppler  orbit  compu¬ 
tation. 

The  precise  orbits  are  computed  from  observations  collected  by  a  worldwide  network  of  Tranet/ 
Magnavox  1502  D.S.  Doppler  tracking  stations  directed  by  DMAHTC  DMAHTC  collects  and 
preprocesses  Doppler  data  and  sends  "clean”  observations  to  N'SWC  The  orbits  are  generated  for  2- 
day  time  spans  by  the  CELEST  orbit  computation  program,  and  stored  as  permanent  files  for  use  with 
the  altimeter  measurements  The  orbit  computation  procedures  and  algorithms  are  described  in 
Reference  2. 


5.0  ALGORITHMS  AND  PROCEDURES 
5.1  PREFILTERING 

Before  they  can  be  reduced  to  best-estimate  along-track  geoid  heights  and  vertical  deflections,  the 
altimeter  heights  must  be  prefiltered.  Prefiltering  entails  the  separation  of  land/ocean  data  and  the 
elimination  of  incorrect  (spurious)  data  that  may  be  caused  by  telemetry  bit  errors,  land  data,  ice 
data,  and/or  environmental  effects.  The  procedure  for  prefiltering  is  outlined  as  follows. 

1.  Computation  of  parameters  required  for  prefiltering,  such  as  geodetic  latitude  and  longitude, 

satellite  rev  number,  longitude  of  the  ascending  node  of  current  rev  and  next  rev,  and  time  of  the 

ascending  node. 

2.  Separation  of  land  and  ocean  data. 

3.  Boundary  value  tests  for  altimeter  height,  SWH,  and  AGC 

4.  Boundary  value  tests  for  SDR  sigma  height,  SWH,  and  AGC 

5.  Missing  data  test. 

6.  Linear  least-squares  fit  to  height,  SWH,  and  AGC. 

7.  Sea  surface  height  boundary  test  (segmenter). 

8.  Vertical  deflection  boundary  test  (after  smoother). 

9.  Identification  of  time-continuous  data  segments 
5  1.1  Geodetic  Latitude  and  Longitude 


See  Section  5  2 
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5  12  Revolution  Number 


Data  associated  with  a  satellite  revolution  around  the  earth  is  identified  by  rev  number, 
beginning  with  the  number  1 ,  increasing  in  increments  of  one  until  the  end  of  the  mission  One  of  the 
inputs  to  the  data  processing  system  is  a  rev  table  or  file  that  is  maintained  and  updated  with  epochs 
of  known  rev  numbers,  longitude  of  ascending  node,  period  of  satellite,  and  shift  of  longitude  of  the 
ascending  node  associated  with  each  epoch  time 


Procedure:  Read  the  rev  file  and  find  the  last  epoch  time  recorded  prior  to  the  time  of  interest  and 
compute  rn,  the  rev  number  of  interest. 


Input: 

re  =  rev  number  at  epoch 


Pe  =  period  of  the  satellite  in  sec  at  epoch 
At  =  difference  between  current  time  and  epoch  (sec) 
Computation:  Calculate  the  new  rev  number,  rn. 


At 

r  ~  r  +  — 
n  e  p 

e 


5  13  Longitude  and  Time  of  the  Ascending  Node 


The  longitude  and  time  of  the  ascending  node  are  computed  using  inputs  from  the  rev  epoch  table 
(Section  5.1.2).  The  algorithms  are  indicated  below. 


Procedure 

Values  in  the  rev  epoch  table  corresponding  to  the  last  epoch  time  recorded  prior  to  the  time  of 
interest  are  used  in  the  computation  of  the  longitude  and  time  of  the  ascending  node 


Input: 


-  rev  number  at  epoch 
=  new  rev  number 

=  shift  in  longitude  of  the  ascending  node,  An  at  epoch  in  deg 
=  An  at  epoch  in  deg 
=  time  of  epoch  in  sec 
=  period  of  epoch  in  sec 
Computation.  Calculate  longitude  of  the  ascending  node,An,  for  rev,  rn 
,\„  =  A  i  ( A  i  /  A.  2 )  A  > 

where 

A  l  -  Anf.  *  AA,,(r„  r,, ) 


AAn 

Ane 

T, 

P„ 
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A2  =  360  deg 

Computation:  Calculate  time  of  the  ascending  node,  T\_ 

T\n  =  Ai  -  ( A1/A2)  A2 
where 

A  i  =  Te  +  (rn-re)Pe 
A2  =  86,400  sec 
5  14  Data  Editing 

Data  editing  is  performed  to  eliminate  spurious  data  points  caused  by  such  factors  as  bit  errors, 
system  problems,  land/ice  contamination,  and  environmental  conditions  Data  editing  consists  of  a 
series  of  tests  and  a  linear  least-squares  fitting  procedure  The  tests  are  based  on  information 
provided  on  the  SDR  and  experience  gained  from  analysis  of  altimeter  data  from  the  GEOS-3  and 
SEASAT  satellites.  The  boundary  values  of  the  tests  are  input  and  were  determined  during  the 
initial  calibration  and  evaluation  phase 

The  altimeter  has  several  modes  of  operations,  shown  on  page  90  of  the  ICD  document  1  Track  1 , 
2,  3,  or  4  are  the  operational  modes  of  interest  for  NSWC’s  altimeter  data  processing  system  Data 
acquired  during  any  other  mode  is  automatically  written  to  the  land/ice  data  file  Data  collected 
while  the  altimeter  is  in  Tracks  1, 2,  3,  or  4  are  processed  as  follows 

1  Separation  of  Land  and  Ocean  Data  Test  See  Section  5  15 

2.  Radar  Altimeter  Height  (H)  Boundary  Test 

Hnnin  <-H(n)>  Umax 

where  H(n)  are  the  measurements  given  at  the  10  pps  rate  in  the  SDR  data  records,  words  5  14, 
and  Hmin  and  Hmai  are  the  minimum  and  maximum  radar  altimeter  heights  expected  The 
minimum  and  maximum  heights  are  input  values  determined  from  orbit  predictions  U pper  and 
lower  limits  for  the  altimeter  heights  are  also  given  in  the  SDR  header  record  If  the  height  is  out  of 
range,  a  flag  is  set  This  flag  is  used  in  the  linear  least  squares  fitting  procedure  for  dubbing  purposes 
and  for  determining  whether  data  belong  on  the  ocean  or  the  land/ice/other  data  file 

3  The  AGC  Boundary  Test 

AGCmm  <  AGC(n)  >  AGCma, 

AGC(n)  is  given  in  the  SDR  data  records,  words  30  39,  at  the  10  pps  rate  and  AGC,,,,,,  and  AGC„1MX 
are  input  values  determined  during  the  validation  and  calibration  phase  If  the  AGC  is  out  of  bounds, 
a  flag  is  set  This  Bag  is  used  during  the  linear  least  squares  fitting  procedure  for  dubbing  purposes 
and  for  determining  if  data  point  is  acceptable 


4  The  SWH  Boundary  Test  * 

SWHmin  <SWH(n)  >SWHma, 


NSWC  TR  86-149 


SWH(n)  is  given  in  the  SDRdata  records,  words  18-27,  at  the  10  pps  data  rate  and  SWHmm  and 
SWHm.»  are  values  determined  during  the  calibration  and  validation  phase.  IfSWH  is  out  of  bounds, 
a  flag  is  set.  This  flag  is  used  during  the  linear  least-squares  fitting  procedure  for  dubbing,  and  for 
determining  whether  the  data  point  is  acceptable 

5.  Sigma  height,  Sigma  SWH,  and  Sigma  AGC  Boundary  Value  Test 

The  sigma  height,  sigma  SWH,  and  sigma  AGC  values  are  tested  in  a  manner  similar  to  height,  AGC, 
and  SWH  boundary  value  tests. 

6.  Mode  Word  Test 

The  SDR  data  record  contains  a  mode  word  that  describes  the  operational  status  of  the  altimeter.  The 
mode  word  consists  of  three  status/mode  sections  of  10  bits  each.  Status/mode  section  #1  is  examined 
for  the  tracking  status.  Data  are  accepted  as  good  ocean  data  if  the  altimeter  is  in  Track  1 , 2,  3,  or  4  . 
otherwise,  data  is  written  to  the  land/ice/other  file.  Flags  in  status/mode  word  section  #3  are  also 
examined  for  tracking  information.  The  following  flags  must  have  the  indicated  values:  gate  index 
equals  1, 2,  3,  4,  or  5;  ACQ,  ACQTC,  attitude,  DHa,  and  LMax  are  equal  to  zero;  and  the  detect  flag  is 
equal  to  1  These  flags  verify  the  quality  of  the  ocean  data. 

7.  Missing  Data  Test 

Information  on  the  SDR  sequenced  by  major  and  minor  frame  count  and  is  written  in  records  that 
contain  10  minor  frames  each  Although  the  data  is  in  sequence,  complete  records  or  data  points 
within  a  record  could  be  missing  Records  are  missing  when  there  is  a  gap  in  the  frame  count;  data 
points  are  missing  when  the  height  equals  zero,  AGC  equals  zero,  or  SWH  equals  zero  and  the  SWH 
"reasonableness  flag”  is  set  SWH  may  have  zero  as  a  valid  value,  therefore,  the  "reasonableness 
flag”  must  be  examined  to  verify  that  the  SWH  value  is  not  valid 

8.  Identification  of  Time-Continuous  Data  Segments 

During  the  data  editing  process,  the  altimeter  data  is  divided  into  time-continuous  segments,  as 
required  by  the  Kalman  smoother  If  any  of  the  following  conditions  is  found  to  exist,  end  of  segment 
and  begin  new  segment  flags  are  set 

°  Data  over  land/ice  encountered 

AGC  has  an  invalid  value  and  land/ice  exist  in  the  surrounding  records 

AGC  has  an  invalid  value  and  data  is  from  an  ice  area  The  bounds  of  the  ice  areas  are  input 
to  the  procedure  since  they  change  with  the  seasons 

’  Invalid  mode  word 

’  Sigma  height,  sigma  SWH,  and  sigma  AGC  from  SDR  are  out  of  bounds 
’  Data  are  missing  between  records 

°  A  combination  of  the  following  conditions  exists  for  a  specified  number  of  consecutive  data 
points 

Invalid  AGC  and  data  is  not  landyice  data 


Invalid  altimeter  height 
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-  Missing  data  within  a  record 
S.  Linear  Least-Squares  Fitting  Procedure 

A  linear  least-squares  Fitting  procedure  eliminates  spurious  data  points. 

Procedure: 

A  straight  line  is  Fitted  to  a  maximum  of  30  consecutive  heights,  and  heights  greater  than  "sigma  x  a 
multiplier”  from  the  straight  line  are  replaced  with  values  from  the  line.  The  process  is  repeated 
until  all  heights  are  within  the  tolerance  or  a  maximum  number  of  iterations  is  reached.  Sigma  is 
computed  for  each  fit  and  used  with  a  multiplier  obtained  from  input.  Values  tagged  prior  to  the  Fit 
are  not  used  in  the  Fit. 

The  SWH  and  AGC  are  processed  in  a  similar  manner 


fj  =  values  to  be  fit  (1  Sj  Smi 
o  =  Sigma  multiplier 
m  =  number  of  values  to  be  fit 
Computation: 

Fit  a  straight  line  to  the  data  points  fj  at  x}  (lSjrSm)as  given  by  the  equations  below  Data  points  are 
given  at  a  constant  interval,  but  some  points  can  be  tagged. 

Y(xj)  =  C0  +  Cixj  (1) 


Atm)  =  y  J 


B(m)  =  X  / 


F  =  m  B(m)—  A  Am) 


m  m  , 

CQ=  B(m)  £  f  -  Atm )  V  jf\  IF 

1  =  1  J  7=1 

ffi  m  \ 

Cj=  -Atm)  f  +  m  y  jf\  IF 


S£»«: 
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Since  Formulas  2  and  3  are  functions  of  the  number  of  data  points,  they  can  be  precomputed.  If  a 
point,  for  example  (x*,  f*)  is  tagged,  then  these  sums  are  modified  as  follows  for  each  tagged  point. 

A(m)  =  A(m)  -  k 

B(m)  =  B(m)  -  k2 

The  sums 

m  m 

Y  f  and  y  jfj 
j= 1  ;= 1 

should  not  include  fk 

Formulas  4, 5,  and  6  are  used  to  compute  Co  and  Ci,  input  to  Formula  1,  yielding  the  fitted  values. 

11.  Sea  Surface  Height  (SSH)  Bounds  Tests 

The  sea  surface  heights  are  globally  known.  In  general,  they  have  a  consistent  range  of  values, 
except  in  areas  near  the  coast  of  India  and  the  Phillipines.  Values  reach  a  minimum  in  the  Indian 
Ocean  and  a  maximum  in  the  Pacific  Ocean  near  the  Phillipines.  Boundary  values  used  in  the 
segmenter  for  eliminating  sea  surface  heights  out  of  range  in  these  areas  are  as  follows.  In  the 
Indian  Ocean,  the  bounds  are  -125m  ^  SSH  •£  125m.  Latitude  ranges  from  -11.5°  to  20°,  while 
longitude  ranges  from  63°  to  90°.  In  the  Pacific  Ocean,  bounds -100m  s  SSH  S  100m.  Latitude 
ranges  from  12°  to  8°,  while  longitude  ranges  from  123°  to  158°.  The  general  global  bounds  are  -80m 
s  SSH  s  80m.  Sea  surface  heights  outside  these  bounds  are  set  to  the  bound  value. 


12.  Vertical  Deflection  Bounds  Test 

After  vertical  deflections  ( VD)  are  computed  in  the  Kalman  smoother,  the  following  test  is  applied. 
-100  arc  secs  VD  s  100  arc  sec 

Vertical  deflections  outside  of  these  bounds  are  set  to  the  bound. 


5.15  Separation  of  Land/Ice  and  Ocean  Data 

The  SDR  is  composed  of  land/ice  and  ocean  data.  One  of  the  requirements  of  the  processing  system  is 
to  separate  the  two  types  of  data.  The  land/ice  data  are  written  to  tape  and  archived.  The  ocean  data 
are  reduced  to  best-estimated  along-track  geoid  heights  and  vertical  deflections  and  then  sent  to 
DMAAC,  NAVOCEANO  and  other  specified  users.  The  procedure  for  separating  the  two  types  of 
data  is  described  below. 

Procedure: 

A  masking  grid  data  file  generated,  by  an  offline  computer  program,  from  the  World2  data  base  of 
coastline  coordinates  consist  of  land/ocean  latitude  and  longitude  boundaries.  Longitude  boundaries 
are  given  at  1-min  intervals  representing  a  total  of  2 1,600  columns.  Latitude  boundaries  within  each 
column  are  ordered  from  zero  min  (South  Pole)  to  10,800  min  (North  Pole).  The  desired  data  column 
is  retrieved  by  longitude.  The  column  is  then  searched  for  the  first  latitude  greater  than  the  desired 
latitude  If  the  entry  number  of  the  latitude  in  the  column  is  odd,  the  location  is  over  land:  an  even 
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entry  number  is  over  water.  Data  values,  representing  min  of  latitude,  are  stored  in  a  FORTRAN  77 
direct  access  file  keyed  on  minutes  of  longitude  and  can  be  randomly  accessed.  The  actual  check  is 
performed  as  follows: 

1.  The  given  minute  of  longitude  is  used  to  access  the  proper  record  of  the  data  Hie. 

2.  These  latitude  values  are  stored  in  a  singularly  dimensioned  array. 

3.  This  array  is  searched  for  the  first  data  value  greater  than  the  given  latitude. 

4.  If  this  array  element  has  an  odd  index,  the  point  is  over  land.  If  the  index  is  even,  the  point  is  over 
water. 

i 

|  Input: 

A  masking  grid  generated  from  the  World2  Data  Base  of  coastline  coordinates  obtained  from 
NAVOCEANO 

i 

i 

|  Output: 

|  A  land  or  ocean  data  point  as  determined  from  the  permanent  data  file  (10,800  X  21,600)  of 

I  land/ocean  latitude  and  longitude  boundaries. 


5.2  SEA  SURFACE  HEIGHT 

The  height  of  the  geoid  above  the  reference  ellipsoid  (See  Figure  3)  and  subsatellite  position 
(geodetic  latitude  and  longitude)  are  computed  from  Equations  1, 2,  and  3  as  given  below 

h  altimetric  height  above  the  geoid  (km) 

r„  -  distance  of  satellite  from  center  of  reference  ellipsoid  (km) 

N  -  height  of  geoid  above  reference  ellipsoid  (m) 

<J>  -  geodetic  latitude  (deg) 

H  geometric  height  above  the  ellipsoid  (km) 

a  -  semimajor  axis  of  reference  ellipsoid  (km) 

e  -  eccentricity  of  the  reference  ellipsoid 

N  =  H  -  h 

where 
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a(l  — e2) 


t  1  \ 
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and  h  is  corrected  for  systematic  and  environmental  effects  as  described  in  Section  5  3. 


SATELLITE 


REFERENCE 

ELLIPSOID 


FIGURE  3.  HEIGHT  OF  GEOID  ABOVE 
REFERENCE  ELLIPSOID 


tan*o=  J=f=t 

V  X  +V 


starting  value  for  4>,  and 

j  ae'tan 4>^ 

^*+l  \/x^  +  y*'  3  Vl  -Ml  -e2)tan2<t>^ 

iterate  until  |  tan  <J>k  + i  -  tan  4>k  I  ^  t. 


where 


i  =  (1  +tan24>i,)  x  10  9,  then 

I  aeltan^/i 

Jt  *  3»  * 


\  =  tan  — , where 0  ^  A  ^  2u 

I  r  * 


4>i  =  <J)ji  + 1  and  A,  are  computed  for  each  observation  lime,  t, 
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I 

a  =  6378.137  km,  e2  =  2f -  P,  and  f  =  1/298  257223563 
5.21  Subsatellite  Velocity 


* 

) 

t 

) 

i 


+ 


'cos  4>  sin  \N 
cos  <J>  sin  \ 
vsin  <f>  , 


|  where  hj  =  0,  since  it  is  the  station  altitude  above  the  reference  ellipsoid. 

> 

i 

|  Routine  NTERP8,  described  in  "The  User's  Guide  for  NTERP8  Routines”  by  Evelyn  C  Durling, 

Code  K14,  may  be  used  to  compute  the  velocity  at  the  given  subsatellite  position  by  differential 
interpolation. 


5.30  CORRECTIONS 


The  altimeter  measurements  are  corrected  for  both  environmental  and  systematic  errors  as 
shown  in  Figure  4.  The  description  of  each  correction  is  given  in  Table  1 . 


Center  Gravity  GEOSAT 
Dh(cg).  Geometric  Distance 
antenna  electrical 


True  range(HC)-H  *  Dh(fm) -Dh(cal) 
-  Dh(inif)  +  Dh(cg) 

•  Dh(SWH/XI) 


Instantaneous  sea  surface  (electrical) 
SWHC.  sea  state  correction 

Instantaneous  sea  surface  (geometric) 
Tide  correction.  DT  (land  and  water) 
Barotropic  (press,  loading)  corr,  D0 

Approximation  to  geoid  height.  GH 


Sea  Surface  Height  (SSH)--h  — HC 

Geoid  Height  ~  oSH  -  Dy  •  Dyy  *  Dq  *  Dj  —  Dg 

Note:  Sea  state  correction  is  combined  with  tilt. 


FIGURE  4.  ENVIRONMENTAL  AND 
SYSTEMATIC  CORRECTIONS 


5.31  Environmental  Corrections 


Environmental  and  system  corrections  applied  to  altimeter  height  measurements  are  shown  in 

Table  1. 
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TABLE  1.  CORRECTIONS  APPLIED  TO  ALTIMETER  HEIGHT 

MEASUREMENTS 


Corrections  Description 

Dh(fm)  Height  bias  resulting  from  FM  crosstalk  with  the  altimeter  return 

signal 


Dh(cal) 


Systematic  height  bias  derived  from  the  analyses  of  calibration  data 


Dh(init) 


Initial  height  measurement  systematic  bias  determined  through 
prelaunch  tests 


Dh(cg)  Height  bias  estimate  caused  by  an  attitude  shift  due  to  change  in 

center  of  gravity  (eg)  as  fuel  is  expended,  and  the  physical  separation 
betv/een  eg  and  the  electrical  center  of  the  radar  altimeter 

Dh(SWH/XI)  Height  bias,  a  function  of  the  off-nadir  spacecraft  (s/c)  orientation 

and  the  SWH 


Dd 

Dw 

Dt 

Db 

Df 


Dry  tropospheric  distance  equivalent 
Wet  tropospheric  distance  equivalent 
Ionospheric  distance  equivalent 

Barotropic  (pressure  loading)  correction  (positive  ( + ),  if  p<  1013  mB) 

Tide  correction  includes  land  and  water.  High  tides  are  assumed 
positive  (  +  ) 


5.311  Tropospheric  Refraction  Correction.  "The  tropospheric  effect  results  in  an  altitude  change  of 
the  order  of  2  5  m  with  the  predominant  portion  of  the  correction  coming  from  the  dry  component  of 
the  atmosphere.  The  maximum  correction  for  the  wet  tropospheric  effect  will  be  less  than  50  cm.”11 
These  limiting  values  are  consistent  with  those  used  in  the  SE  ASAT  data  processing  systems.  Algo¬ 
rithms  for  the  tropospheric  corrections  are  given  as  follows. 


Dry  Tropospheric  Correction 

The  dry  tropospheric  correction  Dp,  is  computed  from  input  pressure  and  latitude. 


(2  277  -  Oil  cos<b  )P 

D  (cm)  =  - - -  ’ 

d  in 


where 

$  =  Latitude  from  location  data 

P=  Pressure  (mB)  If  P  is  unavailable,  use  1013.3. 


I 


NSWC  TR  86-149 


Wet  Tropospheric  Correction 

The  wet  tropospheric  correction  Dw,  is  computed  from  input  surface  temperature  and  partial  pressure 
of  water  vapor. 

Dw(cm)  =  86400.  Pw  /  (273.  +  T / 

where 

Pw  =  Partial  pressure  of  water  vapor  (mB), 

Ts  =  Surface  temperature,  deg  C. 

If  Pw  is  unavailable,  use  12.272. 

IfTg  is  unavailable,  use  10. 

5.312  Barotropic  Correction  The  barotropic  correction,  Db,  is  computed  from  input  atmospheric 
pressure. 

DB(m)  =  009948 (Pa  1013  0) 
where 

Pa  =  Atmospheric  pressure  (mB)  supplied  by  FNOC. 


5.313  Ionospheric  Refraction  Correction.  The  ionospheric  refraction  model  shown  in  Figure  5  was 
developed  by  Dr.  James  Clynch  and  Dr.  Arnold  Tucker  of  the  Applied  Research  Laboratory, 
University  of  Texas,  Austin,  Texas,  and  modified  by  the  late  Dr.  Ralph  Gibson  of  NSWC.  The  model 
uses  telecommunications  predictions  obtained  from  the  National  Oceanic  and  Atmospheric 
Administration,  National  Geophysical  Data  Center,  Boulder,  Colorado. 

The  ionospheric  correction,  Di  =  AR,  given  by  the  equation  below,  is  illustrated  by  Figure  5. 


c  [rs  N{x,y,z)rdr 

AR  "  f2  I  r  y/T^k2' 
e 

where  k  =  re  sin  z 

For  the  GEOSAT  satellite  altimeter  k  =  0,  since  the  range  is  straight  up  and  down,  the  iono¬ 
spheric  correction,  Dj,  is  given  by  the  equation  below.  The  correction  ranges  between  3  and  15  cm 
with  an  accuracy  of  ±  15  percent 


C  [  r  3 

D  (cm)  =  —  N(x,y,z)dr, 

r 

'  € 

where 

r3  =  distance  from  center  of  earth  to  the  satellite  (km) 
re  =  the  radius  of  the  earth  (kml 
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f  =  transmitted  frequency  (MHz) 
c  =  a  known  constant  ( =  8.061389  x  10  6) 

N(x,y,z)  =  electron  density  (el/cm3) 

The  vertical  electron-density  profile  is  shown  in  Figure  6.  The  E  and  FI  layers  are  modelled  by  a 
parabola;  F2  layer  up  to  HmF2  is  modelled  by  a  fourth  degree  curve;  HmF2  to  Hcrit  is  modelled  by  a 
parabola;  and  the  remaining  three  portions  of  F2  are  modelled  by  a  decaying  exponential  with  decay 
constants  Kj_  K2,  and  K3. 


FIGURE  5.  IONOSPHERIC  REFRACTION 
MODEL 

Parameters  for  Ionospheric  Electron-Density  Profile 
F0E  Critical  frequency  of  E  region 

HmE  Height  of  maximum  ionization  of  E  region 

YmE  Semithickness  of  E  region 

F0F1  Critical  frequency  of  FI  region 

HmFl  Height  of  maximum  ionization  of  FI  region 

YmFl  Semithickness  of  FI  region 

F„F2 
H„F2 


Y„,F2 


Critical  frequency  of  F2  region 

Height  of  maximum  ionization  of  F2  region 

Semithickness  of  lower  and  middle  F2  region 
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10.  The  data  smoothing  diagram  is  shown  in  Figure  7.  In  the  filter  or  smoother  model  being  used, 
no  steady  state  is  assumed  for  the  filtering  matrices  This  practice  is  followed  because  of  the 
possibility  of  bridging  over  missing  points.  At  missing  points,  the  filter  procedure  changes  in  the 
error  of  observation,  l/o^2  =  0  This  is  the  add-on  term  in  Dr  Cohen’s  earlier  algorithm  for 
computing  B  from  L  both  forward  and  backward.  At  other  points,  o^2  equals  the  assumed  noise 
variance.  In  the  procedure  in  Section  5  44,  Dr  Cohen  did  not  identify  L  separately,  but  it  is  computed 
in  the  same  way  as  shown  in  earlier  formulation. 


FORWARD  COMPUTATION 
TRANSIENT 
LENGTH 


BRIDGEABLE 


-STEADY  STATE- 


h - STEADY  STATE - 

BACKWARD  COMPUTATION 


TRANSIENT _ I 

'  LENGTH  “I 


GAP  l-_ 

TRANSIENT 

STEADY 

LENGTH 

STATE  - 

L. 

STEADY  _  1  _  ' 

TRANSIENT 

FIGURE  7.  DATA  SMOOTHING  DIAGRAM 

Points  greater  than  a  certain  distance  from  either  of  the  ends  or  from  a  bridged  region  may  be 
assumed  to  be  steady  state  That  is, in  the  forward  direction,  it  could  be  assumed  to  be  steady  state  at 
some  point  from  the  beginning  of  valid  data  continuously  to  the  end  of  valid  data.  At  the  first  point  in 
a  gap,  the  B  matrix  may  be  computed  from  the  previous  one  with  de-weighting,  as  outlined  (Section 
5  44)  by  Dr  Cohen  The  contribution  of  raw  data  to  b  is  also  zero  in  the  gap.  Similar  computations 
are  made  all  the  way  across  the  gap  Likewise,  in  the  backward  direction  at  some  point  from  the  end, 
a  steady  state  may  be  assumed  and  continued  until  the  data  is  interrupted 

The  problem  with  assumption  of  steady  state  is  that  the  distance  required  is  a  function  of  the 
filter  parameters.  Thus  it  is  likely  that  the  worst  case  must  be  assumed  Previously  a  constant 
autocorrelation  distance  permitted  steady  state  to  be  assumed  in  approximately  1 1  sec  The  time 
required  to  arrive  at  steady  state  could  be  an  input 

11.  Splice  with  adjacent  clusters  Adjacent  clusters  are  spliced  together  after  smoothing  with  the 
same  function  as  the  fit  removal  weighting  function.  The  difference  is  that  for  this  splicing  only 
about  20  points  of  overlapping  data  will  be  needed.  That  is  to  say,  20  points  should  be  in  both  clusters 
when  smoothed  As  a  practical  approach,  the  overlapped  20  points  can  all  be  output  with  the  second 
cluster  involved  in  the  overlap 

12.  Restore  fit  removed  and  time  derivative 

13.  Convert  resultant  time  derivative  to  arc  sec  The  X  vector  shown  in  Section  5  44  by  I)r  Cohen  is 
X(  1 )  =  The  geoid  height 

X(2)  =  The  time  derivative  of  geoid  height 
X(3)  =  The  second  time  derivative  (not  needed) 

X(  1 )  is  straightforward  and  ready  for  the  removed  fit  value  to  be  added  to  it 

X(2)  has  the  same  units  as  the  derivative  of  the  removed  fit  curve  (a'2  +  2at;)  +  3  a  t  -  4 )  Const* 
quently,  the  slope  can  be  restored  at  the  point  where  X( 2)  becomes  available  To  obtain  defied  ions  of 
the  vertical  in  arc  sec,  multiply  the  sum  of  the  two  terms  by  206  2648062/velocity  ,  where  velocit  y  is 
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p 

Fi 

(V* 

ftr * 


in  km/sec  Note  that  a  minus  sign  is  present  since  a  positive  slope  is  a  negative  deflection  of  the 
vertical 

14.  Compute  statistics  for  output  header 

15  Write  output 

16  Go  to  2 

5  42  Adaptive  Algorithm  The  computations  for  the  adaptive  filter  parameters,  S\  and  o\  appear  in 
this  section 

N  =  Raw  data  height  (m) 

V  =  Subsatellite  velocity  (km/sec) 

At  =  Time  between  two  Nj  values 
on2  =  Variance  of  geoid  heights 

SN  =  Autocorrection  distance  in  km  Autocorrelation  distance  =  distance  for  autocorrelation 
to  fall  to  1/e  of  value  when  no  separation 

05  =  Expected  noise  of  data 

1  Remove  mean 


> 


i 


ft 


I 


-  . 


N  =  —  ^  N  =  mean  eeoid  height  (m) 
N 

j-  1 


N  -  N  -  N,  l  <  j  <,  N 

j  J 


1  ~ 

L  =  —  ^  A/2,  ACF(O)  •=  LKiriancetm*) 

0  —  j' 

J  =  1 


Set  AX  =  200  km 
2  Let 


K  - 1  NT 


AX 


1 


V  A  7' 
\  h 


A  .S’  =  KVSt 


.  ^  .V  ,V  A  (  T 1  A.S  1  ( m *  1 

i  S-K  — 


1 


\  !A 


\ 


’  N-2K  — 


V  V 


!A 


\(T  1 2 A.S )  im'l 
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[feither  LyL^or  L^L^are  less  than  - ,  let  AX  =  -  AX  and  go  to  2  Continue  calculating  ACF  (iAS),  i  >  2 

i  N~,K  j 

L  =  -  v  N  N  until  L  fl.  <  - 

‘  N-iK  J  J*lK  <  0  e 

7  =  1 

3.  Interpolate  for  Sn 
1 


SN=  1  " 


-L.-L 

p  0  i 

iL  1 _ 

V 

c 

-  AS 

f  J 

y. 

L  .  —  Ld 

i-l  i 

iuS 

4  Compute  On2 


2  90463 


P  (wiener)  = 


X  =  PAS 


<1>(X)  =  (1  +  X  +  -Xle' 
3 


A2  _  ^1 
°N  ~  4><jn 


5  Correct  Sn 


SN=SN  1  +  F(X)  +  0  6117 


(*t)K 


where 


F  (X)  = 


17.74 


/  95  6740  396,000  . 

'  x2  X4  X 6 


A  A  A  2  90463 

X*  =  PCL-SJ.P  =  — x -  l-  =  "AX 

S.v 

n  =  total  number  of  points  in  array  being  processed 
AX  =  VAt 


1  ~  A 

S.v  ~  2  {SN  +  Sv’ 


If  SN  <  80  ,SN  =  80 

A 

•Note  IfX  s5,letS\  =  3Sn 


'  ‘I 
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2  90463 


X  =  PAS 


AiX)  =  (1  *•  X  ♦  ~  X 
3 


2,  -A 


0  852 


o  v  =  1 19  087  o  VP  Wiener  Beta 

N  values  are  expressed  in  meters  Therefore,  variances  computed  from  them  are  in  m2  For  that 
reason,  °s  has  to  be  brought  in  as  meters  However,  when  the  filtering  matrices  are  computed,  all 
three  parameters  must  be  expressed  in  kilometer  units.  S\  will  always  be  expressed  in  kilometers 

5  43  Procedure  for  Removal  of  Cubic  Fit 

1  Generally,  prior  to  computing  adaptive  filter  parameters,  a  cubic  fit  is  removed  from  the  data 
However,  if  there  were  a  short  segment  (less  than  20  points)  not  contiguous  to  other  data,  a  linear  fit 
would  be  sufficient  for  that  case 

2  The  manner  of  removing  the  fit  is  as  follows 

Having  divided  the  data  into  segments  of  approximately  150  sec  each,  take  the  segments  by  pairs  and 
find  the  least  squares  cubic  fit  values  and  slope  values  for  the  pair 

Forj  =  1,1.  where  L  =  length  of  fit  array 

FitZti.j)  =  aj  a 2 1 ( j )  a3t2(j)  +  a4t3(j) 

Slope  W(i,j)  -  a^  +  2a3t(j)  +  3a,|t2(j) 


i  is  the  fit  index 

Next  move  over  one  segment  so  that  the  (i  +  1  )th  segment  is  included  in  both  fits  and  make  a  similar 
calculation 

Use  the  results  of  both  fits  on  the  ( i  +  1  )th  segment  to  find  for  j  =  1 ,  n  where  n  =  segment  length  in 
data  point 

R  <  i  *  1 ,  j  l  "  Z  1 1 ,  n  +  j )  Fiji  *  Z  ( i  1 ,  j  1 1 1  F  (j )  1 

S  ( i  +  1 ,  i )  -  W  ( i ,  n  *  1 1  F  ( i )  *  W  1 1  *  1 .  i )  1 1  -  F  <  i )  I  * 

6  •> 

-  \7.u,  n  *  /)  —  /  ((  *■  \  ,j)  1  (g"  —  gi 

<  n  —  1  i  Ar 
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F'j)  =  1  3g2  -t-  2g3 

At  =  time  interval  between  adjacent  data  points 


J- 1 


3  Subtract  the  fit  values  (R)  from  the  raw  data  values  These  are  the  values  for  processing  from  that 
point  on,  until  after  smoothing,  when  R  and  S  values  are  to  be  added  to  the  smoothing  results 


4  The  double  fit  removal  procedure,  as  shown  in  Figure  8,  assures  continuity  of  the  removed  fit  and 
removed  derivative  Thus,  when  adjacent  segments  are  clustered  and  smoothed,  the  restoration  of 
the  fit  and  derivative  will  not  introduce  a  discontinuity  in  the  smooth  geoid  height  or  deflections  of 
the  vertical  This  procedure  should  also  make  for  a  smoother  transition  when  ends  of  clusters  are 
spliced  together,  when  going  to  a  new  cluster 


ith  fit 


Seg  K 


Seg  K  -  1 


ti  -  1)th  fit - 

Seg  K  -  2 


-t 


FIGURE  8.  OVERLAPS  FOR  CUBIC  FIT  REMOVAL 


5  If  the  data  for  which  the  fit  is  computed  is< ,  the  length  of  one  segment  only,  no  adjustment  to  the 
fit  is  necessary  For  the  first  and  last  segments  of  a  long  span,  no  adjustments  are  made  since  no 
overlapping  fit  exists 

6  The  mean  of  each  segment  is  not  necessarily  zero  after  fit  removal,  since  more  than  one  'egment 
was  involved  Thus  it  will  still  be  necessary  to  remove  the  mean  from  the  segment  for  computing 
filter  parameters 

5  432  Cubic  Fit  Removal  From  Altimetry  Data 
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The  partial  with  respect  to  each  A  must  be  zero 
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V  aAJ2  +  2AT+2AT3  +  2AT*-2TY)  =  0 

2  1  J  -4 

(2A3T*+2AJ2  *  2AJ2+  2  AJb-  2  T2V)  =  0 


Since  T,  =  k,AT.  replace  T  with  k,  to  obtain  A  array  a*  function  of  integers 
Note  IfT  is  made  symmetrical  (center  point  -  01,  then 

S  T  ,  Y  r'.andY  I  s 

lilt  II 

are  all  zero  This  makes  the  <Hi  matrix  easier  to  invert  analytically  and  the  resultant  expressions 
will  involved  fewer  terms 

The  A’s  are  computed  by  removing  powers  of  AT  so  that,  instead  of  summing  over  1  ’  V*  and  l,*' 
the  summing  is  done  over  the  integer  coefficients  to  t he  appropriate  powers 
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Generation  of  Bf,  and  bf 

By  recurrence,  generate  1  Bf,  and  bf)K  and  store  Bf  is  3x3  symmetric  and  bf  is  3x  1 
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Here.Cj  -  I  0  I  and  o«2  is  the  variance  of  the  error  of  observation 

\0/ 

If  the  observation  is  missing 
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Backward  generation  of  solutions 
By  recurrence  backward,  generate  solutions. 

(/Wa2  =  (03’01 

Where  K  =  n>  is  the  last  of  the  observation  times 
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6.0  DATABASE  MANAGEMENT  SYSTEM 


There  are  several  types  of  data  files  involved  in  the  GEOSAT  altimeter  data  processing  system. 
These  files  are  archived,  distributed  to  users,  and  serve  as  inputs  for  geodetic  studies. 

The  magnitude  of  GEOSAT  data  demands  a  systematic  procedure  for  keeping  track  of  it.  After 
several  data  base  management  systems  were  investigated  and  System  2000  was  selected  Subse¬ 
quently,  the  RIM  data  base  management  system  became  available,  and  proved  itself  more  user- 
friendly  as  well  as  needing  less  application  software 

The  basic  yequirement  is  to  maintain  a  log  file  from  which  the  status  of  data  files  can  be  obtained 
In  this  way.  any  specific  information  may  be  retrieved,  written  to  an  output,  and/or  printed 
Provisions  for  adding,  deleting,  and  modifying  entries  are  also  required  The  log  information  consists 
of  a  Revs  Table  and  a  Days  Table,  as  described  in  Section  6  1. 


6.1  LOG  INFORMATION 

The  log  information  appears  as  two  tables,  the  Revs  Table  and  the  Days  Table  One  record  for 
each  satellite  revolution  appears  in  the  Revs  Table  and  one  recored  for  each  day  appears  in  the  Days 
Table.  Tables  2  and  3  list  the  information  contained  on  each 


TABLE  2.  INFORMATION  CONTAINED  BY  THE  REV  TABLE 


Description 

Type 

Format 

Satellite  rev  no 

1 

xxxxx 

Start  year  &  day 

A 

xxxxx 

Start  time  of  rev  in  sec 

R 

\xxxx  XXX 

End  year  and  day 

A 

XXXXX 

End  time  of  rev  isec) 

R 

XXXXX  XXX 

Date  processed 

A 

xxxxx 

Data  quality 

I 

X 

Start  longitude  (deg  ) 

R 

XXX  XXX 

Start  latitude  (deg  ) 

R 

XXX  XXX 

End  longitude  (deg  i 

R 

XXX  XXX 

m 
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TABLE  2.  INFORMATION  CONTAINED  BY  THE  REV  TABLE  (Cont.) 


End  latitude  (deg.) 

Number  of  segments 
Longitude  of  ascending  node 
Year/day  of  ascending  node 
Sec  of  day  of  ascending  node 
Altimeter 
Orbit  adjust 


xxxxx.xxx 


TABLE  3.  INFORMATION  CONTAINED  BY  THE  DAYS  TABLE 


Description 

Year/day 
First  rev  in  day 
Last  rev  in  day 
Original  SDR  tape  ny  # 

Archive  SDR  tape  ny  # 

Date  IGDR  tape  sent  to  AFGL 
AFGL  tape  #  for  IGDR 
Archive  ocean  IGDR  tape  ny  # 
Date  ocean  IGDR  sent  to  NRL 
NRL  tape  #  for  land  IGDR 
Archive  land  IGDR  sent  to  NRL 
Date  land  IGDR  sent  to  NRL 
NRL  tape  #  for  ocean  IGDR 
Archive  FGD  tape  ny# 

Date  FGD  date  sent  to  DM  AAC 
Tape  #  of  FGD  sent  toDMAAC 
Date  FGD  data  sent  to  NAVO 
NAVO  tape  #  for  FGD 
Climatology  data  status  (m.a) 


6.1.1  Queried  Information 


Format 


The  following  operations  may  be  performed  on  the  log  base:  retrieve,  print,  write  to  permanent 
files,  delete,  modify,  and  update. 

Retrieve  and  Print  (and/or  Write  to  a  Permanent  File) 

•Rev  records  for  a  particular  interval 
Rev  records  processed  during  a  particular  time  interval 
■Rev  records  with  longitude  of  ascending  node  within  a  tolerance 


Summary  Given  in  Report  Format 

Total  number  of  revs  ( and/or  list  of  rev  numbers) 
Total  number  of  revs  for  a  part  icular  i  imo  interval 


Total  number  of  revs  sent  to  AC 
Total  number  of  revs  sent  to  N  A  Y< ) 
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Total  number  of  revs  processed  during  a  Total  number  of  SDR  sent  to  AFGL 

particular  time  interval  List  SDRs  by  tape  number  and  revs 

Total  number  of  SDR  by  tape  number  and  contained  on  tape 

revs  contained  on  tape 

Delete,  Modify,  and/or  Update  Entries  in  the  Data  Base 

6.2  GEODETIC  INFORMATION 

"Geodetic  information”  refers  to  the  best-estimated  along-track  geoid  heights  and  vertical 
deflections  contained  on  the  FGD  file.  These  files  are  uniquely  identified  by  rev  number  and  are 
stored  sequentially  on  archival  tapes  Information  regarding  these  files,  such  as  their  location  and 
status  can  be  extracted  from  the  log,  computer  programs  manipulate  this  data  as  required 

In  addition  to  the  distribution  of  complete  rev  files  to  users,  there  will  be  some  requests  for  rev 
data  located  in  a  particular  area  of  the  world.  This  requires  specialized  data  base  management 
procedures.  In  the  future,  this  problem  will  be  treated  in  greater  detail. 

7.0  UTILITY  SOFTWARE 

The  major  software  for  this  system  generates  along-track  geoid  heights  and  vertical  deflections. 
Additional  software  permits  such  peripheral  type  processing  as 

writing  user  compatible  tapes  archiving 

•reading  tapes  generated  by  other  computers  special  plotting 

converting GEOSAT  files  for  use  in  existing  programs  determining  revs  in  an  area 
extracting  rev  data  in  an  area  validating  data 

8.0  GEOSAT  DATA  FILES,  SOURCES,  AND  USERS 


Data  File  Name 


Source  (Computer 


User  (Computer) 


Sensor  Data  Record 

JHU/APUSEL  32/77) 

NSWC  (C DC/865 
and  875) 

Average  Data  Record 

JHU/APUSEL  32/77) 

NORDA 

Intermediate  Geophysical 

Data  Record 

NSWC  (C DC/865  and  875) 

NRL  (Data  General 
MV80000) 

SDR/Precise  Orbit 

NSWC  (CDC/865und  875) 

AFGL  (CYBER  72 
CDC6600  VAX  11  780) 

Unclassified  Data  Record 

Civilian  Community 

Filtered  Geophysical  Data 

NSWC  (CDC/865und  875) 

DMAAC 
(UNI  VAC  1 162) 

& 


I 


Vlt* 
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Data  File  Name 


Source  (Computer 


User  (Computer 


Filtered  Geophysical  Data 


N'SWC  (CDC/865and  875) 


NAVOCEANO 
(VAX  or  L'NIVAC  1182) 


Environmental  Data 


FNOC  (CDC/NODDS) 


NSWC  (NODDS/Mohawk) 


9.0  NEW  GENERAL  PURPOSE  COMPUTERS 


During  development  of  the  GEOSAT  analysis  and  data  reduction  system,  NSWC  acquired  new 
general  purpose  computers,  C DC/865  and  875.  Unclassified  operations  are  performed  on  the 
CDC/865;  classified  operations  on  the  CDC/875.  While  the  computers  have  ample  storage  space  for 
handling  the  huge  number  of  GEOSAT  data  files,  cost  prohibits  their  use  for  archival  purposes.  Files 
are  copied  to  tapes  for  archival 


10.00  DATAFILE  FORMAT  AND  CONTENT 


File  formats  and  contents  for  the  SDR,  1GDR,  FGD  are  listed  in  this  section,  as  well  as  precise 
Doppler  orbit  data. 

10.1  SDR 

The  SDR  tapei  is  among  the  data  products  generated  at  APL  It  consists  of  altimeter  data  for  a  1  - 
day  (24  hr)  period  and  has  one  header  record,  followed  by  a  series  of  data  records  given  at  1-sec 
intervals. 


10.2  IGDR 

The  IGDR  is  an  intermediate  file  generated  by  the  preprocessor  section  of  the  data  processing 
system.  The  file  is  basically  the  SDR  with  the  addition  of  such  information  as  the  satellite  rev 
number,  height  above  the  reference  ellipsoid,  longitude  and  time  of  the  ascending  node,  geodetic 
latitude  and  longitude,  and  environmental  data.  The  file  consists  of  the  header  record  taken  from  the 
SDR  file,  rev  header  record,  rev  data  records  for  each  second  of  data,  and  an  end  of  file.  The  previous 
two  record  types  are  repeated  for  all  revs  processed  from  the  SDR.  The  format  and  contents  of  the 
different  records  are  given  below 

10.21  IGDR  Header  Record  Contents 


Format 


Approximate  Range 
of  Significance 


Description 


First  word 


SDRxxvvdddnn 


SDR  Tape 


SDR  HEADER 


Identification 


NSWC  TR  86-149 


j  Word 

s 

Format 

Approximate  Range 
of  Significance 

DescriDtion 

Last  word  of  SDR  +  1 

A8 

IGDRxxnn 

Identification  for  IGDR  file 
xx  =  Code  Version 
nn  =  Tape  Version 

i  "2 

15 

yyddd 

Processing  date  of  IGDR 

F7.6 

xxxxxx 

Time  interval  (sec) 

>.  + 4 

* 

St 

\ 

\ 

k 

F8.3 

xxxx.xxx 

Semimajor  axis  axis  WGS  84 
Ellip  (km)  for  use  in 
data  reduction 

j 

F7.3 

XXX. XXX 

Flattening  (reciprocal) 

WGS84 

+  6 

F4.3 

XXX 

Timing  correction  (sec) 

S  +7 

E 

11 

X 

0  =  ocean  data 

'  +7 

• 

1  =  land  data 

f 

i  10  22  IGDR  Rev  Header 

\ 

|  Word 

f 

Format 

Approximate  Range 
of  Significance 

Description 

J  1 

> 

15 

XXXXX 

rev  number 

/ 

<  2 
< 

12 

y'y 

Year  (first  observation  of  rev) 

1  3 

W 

13 

ddd 

Day  (first  observation  of  rev' 

W 

4. 

FI  1  6 

XXXXX. XXXXXX 

Sec  of  day  (first  observation 
of  rev) 

t 

F6.5 

X. XXXXX 

Longitude  of  the  ascending 
node(0  to  2  n,  EOG)  (RD) 

6. 

12 

yy 

Year  of  longitude  of  ascending 
node 

!  7. 

13 

ddd 

Day  of  longitude  of  ascending 
node 

8 

FI  1  6 

XXXXX. XXXXXX 

Time  of  equator  crossing 
(sec  of  day) 

)  9 

f 

j 

F6.5 

X  XXXXX 

Longitude  of  the  ascending 
node  of  next  rev 
(0  to  2  n,  EOG)  (RD) 

• 

• 

40 
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Approximate  Range 

Word 

Format 

of  Sienificanee 

DescriDtion 

10. 

12 

yy 

Year  of  longitude  of  ascending 
node  of  next  rev 

11. 

13 

ddd 

Day  of  longitude  of  ascending 
node  of  next  rev 

12. 

F 1 1  6 

xxxxx.xxxxxx 

Time  of  equator  crossing 
of  next  rev 

13. 

18 

xxxxxxxx 

Major  frame  count  from  SDR 
(first  observation  of  rev) 

14. 

12 

XX 

Minor  frame  count  from  SDR 
(first  observation  of  rev) 

15.  F7.3 

10.23  IGDR  Data  Record 

XXXX. XXX 

Approximate  Range 

Period  of  satellite  (sec) 

Word 

Format 

of  Significance 

Description 

1 

15 

xxxxx 

rev  number 

2. 

12 

yy 

Year  of  first 

observation  in  record 

3 

13 

ddd 

Day  of  first 

observation  in  record 

4. 

Fll  6 

xxxxx.xxxxxx 

Sec  of  first 

observation  in  record 

5. 

no 

xxxxxxxxxx 

Mode  from  SDR  Data  Record 

6. 

110 

xxxxxxxxxx 

•FLAGWORD  1  (AGC, 
Height, Land) 

7. 

110 

xxxxxxxxxx 

•FLAGWORD  2  (Mode, 
Missing  Data,  rev) 

8 

F9  6 

xxx.xxxxxx 

Height  ( 1 )  from 

Altimeter  (KM)) 

17 

F9  6 

XXX  xxxxxx 

Height  (10)  from 

Altimeter  (KM) 

18 


14 


xxxx 


Sigma  H  (MM) 
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Format 


Approximate  Range 
of  Significance 


Description 


r ©s. 


±  xxxx 

Delta  H  (SWH/off-Nadir  ang) 
(MM) 

^  < 

±  XXX 

Delta  H-FM  crosstalk  (MM) 

±  x.xxxxx 

Latitude  ( 1 )  ( Rad)  (  ±  rx/2) 

±  x.xxxxx 

Latitude  (10)  (Rad)  ( ±  n/2) 

x  xxxxx 

Longitude  (1 )  (Rad)  (0-to  2, 
EOG) 

x.xxxxx 

Longitude  (10)  (Rad)  (0  to  2n, 
EOG) 

4 

xxx.xxxxxx 

Height  (10)  above  reference 
ellipsoid  (KM) 

-v 
■> 
■r  w, 

xxx.xxxxxx 

Height  ( 1 0)  above  reference 
ellipsoid  (KM) 

•  ."j 

r  -  - 

• 

t  • 

±  XX. XX 

SWH  (1)  (M) 

« 

±  X  XX 

SWH  (10)  (M) 

Ow" 

vl 

.*•  .*■ 

V', 

X  XX 

♦Sigma  SWH  (M) 

V- 

•V.-. 

±  X  XX 

Delta  SWH 

(SWH/ofT  Nadir  Ang)(  M ) 

C 

±  XXX. XX 

AGC(IHDB) 

vN 

V.*J 

±  xxx  xx 

AGC  (10)  (DID 

i 

-  4 

csvi 

XX  XX 

tx  XX 

♦Sigma  AGC  (DID 

Delta  AGC  (SWH  off 

Nadir  Ang)  ( DB) 

* 

i  x  xx 

Delta  AGC  HIM DB) 

-  V 

t  X  XX 

Delta  AGC  iTl(I)B) 
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Approximate  Range 

Format  of  Significance  Description 


77.  F3.0  ±  xx.  DH/DTIM/S) 

78  F3.2  x  xx  Off-nadir  angle  ( Deg) 

79  F4.2  xx  xx  Sigma  zero  (backscatter 

coefficient)  (dB) 

80  F3.1  xx  x  Windspeed  (W  (C))  (M/S) 

81  F4.3  xxxx  Vatt-subT(V) 

82  F4.1  xxx .  x  RCVTD  (Deg  C) 

83  P3 . 1  xx  x  Ionospheric  correction 

(CM) 

84  F3.1  xx  x  Wet  Tropospheric 

Correction(CM) 

85  F4.1  xxx  x  Dry  Tropospheric  Corr  (CM) 

86.  F4.2  ±x.xx  Barotropic  Correction  (M) 

87  F5.2  ±xx  xx  Tide  Correction  (M) 

88  14  xxxx  Atmospheric  Pressure 

(MB) 

89  F4.3  x  xxx  Average  Velocity  (K M/S) 

*  From  NSWC  Preprocessor  program  (30-point  fit) 

10.3  FGI)  RKCORI) 

The  FGD  file  is  the  final  output  of  the  NSWC  altimeter  data  processing  system.  The  file  contains 
information  for  one  revolution  of  the  satellite  and  is  divided  into  time-continuous  segments  of  data 
based  upon  autocorrection  distance  Kach  segment  is  composed  of  a  header  record  and  duta  records 
for  each  data  point  in  the  segment  Segments  are  separated  by  an  EOF 

10  31  FGD  Header  Record  Contents 


Format 


Approximate  Range 
of  Significance 


Description 


XXXX 

number  of  points  in  segment 

xxxxx 

rev  number 

vyddd 

year/date  data  processed 

V.V-  /xV>_ 


NSWC  TR  86-149 


Format 


Approximate  Range 
of  Significance 


Description 


XX 

starting  year  of  segment 

XXX 

starting  day  of  segment 

XXXXX. XXX 

starting  sec  of  segment 

xxxxxx 

time  interval  (sec) 

xxxx. 

autocorrelation  distance  (km) 

XX. XX 

standard  deviation  of  data 

(m) 

XX.  XX 

standard  deviation  of 
geoid  heights  (m) 

XX. XX 

standard  deviation  of  vertical 
deflections  (arc  sec) 

x.xx 

rms  of  filtered-raw 

differences  (m) 

xxxx 

average  subsatellite  velocity 
(km/s) 

±  xxxx 

height  calibration  bias  (mm) 

±  xxxxx 

system  height  bias  (mm) 

xxxxx 

height  bias  due  to  change  in 
center  of  gravity  (mm) 

XXX 

time  correction  (sec) 

XXXXXX  XXX 

speed  of  light  (km/sec) 

XXXX  XXX 

semimajor  axis  WGS  84  ref. 
ellip  (km) 

XXX. XXX 

flattening  reciprocal-WGS  84 

XXXX. XXX 

satellite  period  (sec) 

X  XX 

orbit  uncertainty  in  geoid 
heights  (m) 

XXX 

orbit  uncertainty  in  vertical 
deflections  (arc  sec) 

X  XXXXX 

longitude  of  ascending 
node (rad) 

w  Wt 
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Word 

Format 

Approximate  Range 
of  Siznificance 

DescriDtion 

25 

12 

XX 

year  of  longitude  of  ascending 

26 

13 

XXX 

node 

day  of  longitude  of  ascending 

27 

F9  3 

XXXXX  XXX 

node 

sec  of  day  of  longitude  of 

28 

F7  5 

X xxxxx 

ascending  node 

longitude  of  ascending  node 

29 

12 

XX 

of  next  rev  iradl 

year  of  longitude  of 

30 

13 

XXX 

ascending  node  of  next  rev 

day  of  longitude  of  ascending 

31 

F9  3 

XXXXX. XXX 

node  of  next  rev 

sec  of  day  of  longitude  of 

A  file  is  the  consolidation  of  groups  of  150  sec  (or  lk/km)  of  time 

ascending  node  of  next  rev 

-continuous  filtered  data  with 

approximately  the  same  autocorrelation. 

10  32  FGD  Data  Record  Contents 

Word  Format 

Approximate  Range 
of  Sienificance 

Description 

1 

F9  3 

XXXXX. XXX 

time  (sec) 

2. 

F8  5 

±  X  XXXXX 

latitude  (rad) 

3 

F7  5 

X  XXXXX 

longitude  (rad) 

4 

F7  2 

±  XXX  XX 

geoid  height  (m) 

5 

F6  3 

xxxxx 

•sigma  h  (m) 

6 

F7  2 

±  xxx  xx 

vertical  deflection  (arc  sec) 

7. 

F7  2 

±  xxx  xx 

raw  geoid  height  uncorrected 

8 

F7  2 

±  xxx  xx 

(m) 

raw  geoid  height  with 

9 

15 

xxxxx 

correction  applied  (m) 

flagword  (see  table) 

±  XX  XX 


significant  wave  height  (mi 
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Approximate  Range 


Word 

Format 

of  Significance 

DescriDtion 

11 

F4  2 

X  XX 

•sigma  SWH  (ml 

12. 

F7  2 

±  XXX  XX 

AGC  (db) 

13 

F5  2 

XX. XX 

•sigma  AGC  (db) 

Heieht  Corrections 

14. 

F5  0 

±  XXX 

FM  crosstalk  bias  (mm) 

15 

F6 .0 

±  XXXX 

height  (off-nadir,  SWH 
component)  (mm) 

16 

F5  1 

±  XXX 

ionospheric  correction  (cm) 

17 

F4  1 

XX  X 

wet  tropospheric  correction 
(cm) 

18 

F5  1 

XXX  X 

dry  tropospheric  correction 
(cm) 

19 

F5  2 

±  X  XX 

barotropic  correction  (m) 

20 

F6  2 

±  xx  xx 

tidal  correction  (m) 

21 

F4  2 

X  XX 

tilt  (off-nadir  angle)  (dg) 

22 

F4.1 

XX  X 

windspeed  (nv’s) 

23 

F5  2 

XX. XX 

ocean  backscatter  coefT  (db) 

24 

F6  2 

±  XX  XX 

SWH  correction 
(SWH/off-nadir  angle)  (m) 

25 

F6  2 

±  XX  XX 

AGC  correction  (SWH  off 
nadir  angle)  (db) 

26 

F6  1 

±  XXX  X 

steric  anomaly  (cm) 

‘parameters  from  NSWC  preprocessor  program 

The  fiagword  is  a  5-digit  integer  revealing  the  status  of  each  of  the  parameters  listed  below  At 
present,  the  rightmost  bits  of  integer  are  reserved  for  this  purpose  N  umbering  is  from  right  to  left 

Parameter  Status 

1.  Unfiltered  Geoid  Height  l  out  of  bounds  Use 

-125m  S  gh  <  125  m  for  area  (deg) 
lat  1 1  5  long  64  0 
lat  20  0  long  89  0 
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Parameter 


2  Vertical  Deflection 


3  Significant  Wave  Height 


Parameter 
4  Automatic  Gain  Control 


5  Ionospheric 


6  Wet  Tropospheric 


7  Dry  Tropospheric 


8  Harotropic 


9  Tide 


10  Dubbed  In  Data  Point 


Status 


•100m  ^  gh  ts  100m  for  area  (deg) 
lat  12  0  long  123  0 
lat  8  0  long  158  0 

80m  s  gh  ^  80m  for  remaining  area 
0- within  bounds 

1  out  of  bounds  L'se 

(  100  arc  sec  Svd  s  100  arc  sec) 

0- within  bounds 

l  -out  of  bounds  Use 

(0m  <;  SWH  <;  20m) 

0  within  bounds 

Status 

1  out  of  bounds  Use 

<22db  <  AGC  <;  38  db) 

0- within  bounds 

Environmental  Corrections 


1  out  of bounds 

(0cm  <  iono  <  25cm) 

0  model  value 

1  out  of  bounds 

(0cm  :S  wtrop  ■£  44cm) 
standard  value  used  (13  2  cm) 

0  model  value 

1  out  of  bounds 

(194cm  5  dtrop  <  247cm) 
no  data,  use  standard  pressure  value  (1013  3) 
0  model  value 

I  out  of bounds 

(70  4  cm  <  baro  <  1 56  2  cm  i 
no  data,  use  (0cm) 

0  model  value 

1  out  of  bounds 

(  10m  ?  --  tide  <  ^  1  dm) 
no  data,  use  lOcnu 
0  Schwiderski  model 

I  data  dubbed  in  to  fill  a  small  gup  of  missing  data 
0  no  dubbing 


rV'> 


11  Tilt 


1  lilt  value  greater  than  allowed  field  width 
0  within  allowed  field  width 
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10.4  DOPPLER  PRECISE  ORBIT 


The  Cclest  computer  program^  takes  as  input  raw  data  that  is  obtained  by  the  network  of  Dopplc 
tracking  stations  and  generates  a  precise  orbit  Ultimately,  precise  orbits  will  be  generated  for  the 
entire  mission  and  employed  with  altimeter  measurements  to  compute  raw  geoid  heights 

10  41  Propagated  Trajectory  Kile  (earth-fixed  or  inertial)  (Tape  10) 

The  propagated  trajectory  file'^  may  consist  of  one  or  several  trajectories,  depending  upon  the 
running  option  selected.  The  file  contains  two  record  types  Type  1,  a  header  for  each  trajectory  that 
appears  only  once  for  each  trajectory,  and  record  Type  2,  containing  the  position  and  velocity  of  the 
satellite  for  a  given  time.  Record  Type  2  is  repeated  for  each  time  line  for  each  trajectory  The  last 
record  of  a  trajectory  is  a  record  of  Type  2  containing  a  negative  value  for  the  consecutive  record 
number 


FORMAT  OF  PROPAGATED  TRAJECTORY  FILE 
RECORD TYPE  1 


rd  No. 

Type 

Symbolic  ? 

S'ame  Description 

1 

1 

NB 

Record  no  =  l 

2 

F 

Trat(l) 

Year  of  epoch  of  trajectory 

3 

F 

Trat (21 

Day  of  epoch  of  trajectory 

4 

F 

Trat  (3) 

Sec  of  epoch  of  trajectory 

5. 

F 

Trat  (4) 

Tve  to,  time  of  vernal  equinox  minus 
epoch  of  the  trajectory 

6 

F 

Trat  (5) 

Interval  at  which  trajectory  is  written 
(sec) 

7 

F 

Trat  (6) 

hast  time  on  trajectory  (sec  from  epoch) 

8 

F 

Trat  (7) 

Integration  interval  used  when  creating 
the  trajectory 

9 

I 

(flow 

Kind  of  trajectory  =  10,  11,  12,  or  13 

10 

I 

11X1) 

Satellite  number 

1 1 

1 

KDS<1) 

Indicates  if  trajectory  was  made  using 
epoch  of  date  -  1,  or  epoch  of 

1950  0 

12 

I 

IIX3) 

Improvement  cycle  number 

13 


A 


Idate 


Time  dock  value  when  trajectory  was 
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Word  No 

Type 

Symbolic  Name 

Description 

14 

1 

Nrev 

Revolution  number  of  satellite 

15 

1 

None 

Dummy  word 

RECORD TYPE  2 

Word  No. 

Type 

Symbolic  Name 

Description 

1 

l 

NB 

Consecutive  Record  No  =  2,3,  ,  N 

2 

F 

TI 

Ti-To  (sec  from  epoch) 

3 

F 

Tra  (2) 

x  -  earth-fixed  or  inertial  components  of 
satellite  position  at  time  (Ti) 

4 

F 

Tra  (3) 

y  -  earth-fixed  or  inertial  components  of 
satellite  position  at  time  (Ti) 

5 

F 

Tra  (4) 

z  -  earth-fixed  or  inertial  components  of 
satellite  position  at  time  (Ti) 

6 

F 

Tra  (5) 

xdot  -  earth  fixed  or  inertial  component 
of  satellite  velocity  at  time  (Ti) 

7 

F 

Tra  (6) 

ydot  earth  fixed  or  inertial  component 
of  satellite  velocity  at  time  (Ti) 

8 

F 

Tra  (7) 

zdot  -  earth-fixed  or  inertial  component 
of  satellite  velocity  at  time  (Ti) 

10.5  TIME  CORRECTION  TABLE 

Each  SDR  spans  approximately  one  day’s  worth  of  data  Time  tag  information  (year,  day,  sec)  and 
correlated  frame  count  are  provided  in  the  header  record  for  the  beginning  and  end  of  each  data  span 
This  information  allows  accurate  time  to  be  assigned  to  each  data  record  Nevertheless,  analyses 
have  occasionally  turned  up  errors  in  the  SDR  header,  indicating  that  compensatory  time  corrections 
are  required.  Time  corrections  for  the  appropriate  days  are  placed  in  a  table  and  used  as  input  to  the 
altimeter  data  processing  system.  The  format  of  the  time  corrections  data  table  is  given  below 

Year 

Beein  Time  Correction  End  Time  Correction 

.V  .-.  -N  . N  . V’ 
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10.6  SUN  SPOT  TABLE 


The  ionospheric  refraction  model  described  in  Section  5  313  includes  the  sun  spot  number  in 
the  computation  of  the  ionospheric  refraction  correction  to  the  altimeter  height  These  numbers  are 
given  on  a  monthly  basis  and  stored  in  a  table  that  is  input  into  the  altimeter  data  processing  system. 
The  format  of  the  sun  spot  table  is  given  below 

Year  Month  Sun  Spot  Number 

15  F10  3 

10.7  SATELLITE  REVOLUTION  EPOCH  TABLE  TABLE 

The  satellite  rev  number,  longitude  and  time  of  the  ascending  node  are  computed  by  the  prepro¬ 
cessor  program.  The  algorithms  for  computing  these  values  require  accurate  reference  values.  Due 
to  GEOSAT’s  orbital  eccentricity,  the  reference  values  are  updated  every  100  revs  and  stored  along 
with  other  pertinent  information  in  a  table  Included  as  pertinent  information  are  the  satellite  period 
and  the  shift  in  longitude  of  the  ascending  node,  computed  from  the  precise  orbit  and  averaged  over  a 
two-day  period.  The  table  is  input  to  the  preprocessor.  Its  format  is  given  below. 


@111 

Year 

Day 

Longitude  of  Node 

Longitude  Shift 

15 

F12.6  F83 

F9  5 

F9.5 

11.0  ACRONYMSAND  ABBREVIATIONS 


11.1  AGENCIES 

AFGL 

Air  Force  Geophysical  Laboratory 

ARL 

Applied  Research  Laboratory,  University  of  Texas 

CIA 

Central  Intelligence  Agency- 

DMA 

Defense  Mapping  Agency- 

DM  A  AC 

Defense  Mapping  Agency  Aerospace  Center 

DMAHTC 

Defense  Mapping  Agency  Hydrographic  Topographic  Center 

DON 

Department  of  the  Navy 

FNOC 

Fleet  Numerical  Oceanography  Center 

JHU/APL 

Johns  Hopkins  University  Applied  Physics  Laboratory 

NASA 

National  Aeronautical  Space  Administration 

NAVOCEANO 

Naval  Oceanographic  Office 

NO  A  A 

National  Oceanic  and  Atmospheric  Administration 

NORDA 

Naval  Ocean  Research  and  Development  Activity- 

NRL 

Naval  Research  Laboratory 

NSWC 

Naval  SurfaceWeapons  Center 

I 
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11.2  OTHER 

N'DR  NORDA  data  record 

FGD  Filtered  Geophysical  Data 

GOAP  GEOSAT  Ocean  Applications  Program 

IGDR  Intermediate  Geophysical  Data  Record 

pps  Points  per  sec 

rev  One  revolution  of  the  satellite  commencing  at  the  equator  (ascending  node)  and 

ending  at  the  equator 
s/c  spacecraft 

SDR  Sensor  Data  Record 

UDR  Unclassified  Data  Record 

WDR  Waveform  Data  Record 
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